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Abstract 

Consider on a manifold the solution X of a stochastic differential equation 
driven by a Levy process without Brownian part. Sufficient conditions for 
the smoothness of the law of Xt are given, with particular emphasis on non 
compact manifolds. The result is deduced from the case of affine spaces by 
means of a localisation technique. The particular cases of Lie groups and 
homogeneous spaces are discussed. 



1 Introduction 

Consider a M. m - valued Levy process At without Brownian part, a d-dimensional man- 
ifold M, and the M-valued solution X t of an equation 

X t+dt = a{X t , dA t ) + b(X t )dt, X = x , (1) 

for coefficients a and b such that a(x, 0) = x. The precise meaning of this equation 
will be given later, as well as conditions implying the existence and uniqueness of a 
solution. The aim of this article is to give sufficient conditions ensuring the smooth- 
ness of the law of X t at any time t > 0. We also study more precisely the case 
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of Levy processes on Lie groups, and of some classes of processes on homogeneous 
spaces. 

Proving the smoothness of the law of a random variable has motivated, in the case 
of continuous diffusions, the introduction of Malliavin's calculus. When M = W d , 
Bismut's approach to this calculus has proved to be useful for processes with jumps 
which are solutions of equations of type (EQ), and this topic has been intensively 
studied since [5] and [4]. Different techniques, each of them having its own domain 
of applicability, have been introduced afterwards. These techniques can be roughly 
divided into two classes. 

The first class relies, as in [5j, on some infinitesimal perturbations (in space or 
in time) on the jumps of the Levy process A. A differential calculus can be based 
on these perturbations, and the associated integration by parts formula enables to 
study the smoothness of the law of X t . 

The second class of techniques (also when M = M. d ) has been worked out in [24J ; 
rather than a differential calculus, one uses a finite difference calculus consisting in 
appending and removing jumps. This is not a differential calculus so there is no inte- 
gration by parts formula in the usual sense, but there is still a duality formula which 
can be written on the Poisson space of jumps, and this formula can be interpreted 
as a duality between appending and removing jumps, see [221 [23]; this calculus has 
been applied to the smoothness of the law of X t by [24J. Its advantage is that no 
smoothness is required for the Levy measure of A; in particular this measure may 
have a countable support. The proof of the smoothness of the law is based on an 
estimation of the characteristic function (Fourier transform) of X t ; if this function is 
proved to decrease rapidly at infinity, then X t has a C°° density. In [24] (as it will 
be the case in this article), it is assumed that A has no Brownian part and that the 
coefficient a satisfies a non degeneracy condition similar to the ellipticity condition 
for continuous diffusions; notice that the class of so-called canonical equations with a 
Brownian part was studied with this finite difference calculus by [14] , and Hormander 
type conditions were considered by [T7], also for canonical equations. 

The extension of [24] to the case where M is a manifold is not immediate. The 
basic tool of the case M = M. d was the Fourier transform. When M is a symmetric 
space, the Fourier transform can be replaced by the so-called spherical transform, 
and this has been used by [20] in order to prove the smoothness of the law for a 
class of processes (but the Levy measure was not allowed to be purely atomic in the 
case without Brownian part). However, adapting the Fourier transform technique 
of M. d seems difficult for general manifolds. When M is diffeomorphic to M. d , we 
can of course apply the result for M. d , but the assumptions needed for this result 
are generally not translated into canonical assumptions on M: they depend on the 
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choice of the diffeomorphism, and they may be hard to verify, see in Section 15.11 the 
example of the hyperbolic space. Actually, even the case where M is an open subset 
of M. d is not trivial. 

In order to get around these difficulties, our aim here is to apply some localisation 
techniques; we consider an atlas of M and use the results of the affine case on each 
local chart. However, it is known that localisation is made difficult by the presence 
of jumps. In [26], we have applied such a technique in order to study the smoothness 
of harmonic functions on some domain D of M d , and it appears that these functions 
are not always C°°, even for processes with a C°° density; their order of regularity 
depends on the number of jumps needed to exit D, so that the smaller the jumps are, 
the smoother the function is. Our plan is therefore, first, to apply the localisation 
technique and obtain the C regularity of the density when big jumps are removed, 
and then check that adding big jumps does not destroy this smoothness. In order 
to conclude, when the manifold is not compact, we have to make an assumption on 
the size of these big jumps, namely that for any relatively compact subset U of M, 
the set of points from which the process can jump into U is relatively compact (one 
cannot jump from a very distant location). This condition can be viewed as dual to 
the condition of for the smoothness of harmonic functions. The result (Theorem 
[2]) is given in Section [2j after the particular case of canonical equations (Theorem [1]); 
it is proved in Section [31 

Then we relax this "big jumps" condition in three cases: 

• when the Levy kernel for these big jumps is smooth (Theorem [3]), 

• when M is a Lie group G and X t is a Levy process on it, 

• when M is an homogeneous space G/H and X t is obtained by projecting on 
M a Levy process on G. 

Lie groups and homogeneous spaces are the purpose of Section 0J 

In Section [5l we give some examples, and also some counterexamples where the 
"big jumps" condition is not satisfied, and the smoothness of the law fails. 

2 The smoothness result on manifolds 

In this section we give a precise meaning for Equation (TjQ), but before that we intro- 
duce the manifold M and the Levy process A t . Then we give the assumptions on the 
equation, give the smoothness results about the law of the solution, and make some 
comments on our model. Actually, before explaining the general case, we state the 
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main result in the particular case of so-called canonical equations; assumptions are 
indeed much easier to write in this case. The proofs are postponed to Section [31 

2.1 The manifold 

The manifold M is supposed to be Hausdorff, separable, paracompact, C°° and of 
dimension d; we do not suppose that it is connected; it may have an at most countable 
number of connected components, and our processes will be allowed to jump from a 
component to another. The tangent bundle is denoted by TM = [J T X M. If M is 
not compact, we consider its one-point compactification MU {oo}; if M is compact, 
oo is a point which is disconnected from M. This additional point will be viewed as 
a cemetery point; this means that real functions / on M are extended by putting 
/(oo) = 0, and that a process on MU {00} which hits 00 or tends to 00 stays at that 
point forever. 

Under these assumptions, one can embed M into an affine space 1SL N by means of 
Whitney's theorem, and one can consider Riemannian metrics on M by using a par- 
tition of unity; Riemannian distances are of course only defined between points in the 
same connected component, and they otherwise take the value +00. We would like 
to have intrinsic assumptions, which do not depend on a particular embedding or a 
particular metric. This will be possible because Riemannian distances are equivalent 
on compact subsets, so a metric condition which is supposed to hold uniformly on 
compact subsets will actually not depend on the choice of the metric. For instance, 
if / is a C l function and K is a compact subset of M, if we denote by D^f(x) the 
iterated differential of / (this is a multi-linear form on (T X MY), we can define 

1 

\f\ K/ = Vsup|£>V(z)| 

if we have chosen a Riemannian metric on M, but changing the metric will lead to 
an equivalent semi-norm; by allowing K to vary, we obtain a Frechet space C e (M). 
A signed measure v on M is said to be absolutely continuous, respectively C e , if its 
restrictions to local charts are absolutely continuous, respectively have C e densities 
with respect to the Lebesgue measure; it is said to have positive density if it has 
(strictly) positive density with respect to the Lebesgue measure on charts; this does 
not depend on the atlas; let M. (M) be the set of C measures. If we choose a C°° 
reference measure dx with positive density, we can define the family of semi-norms 

\v\k,£ = \dvjdx\ 
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Changing the reference measure and the Riemannian metric will not change the 
topology of M e (M). We also let 

Mi = [u G M e (M); u(K c ) = 0} 

if K is a compact subset of M. 

The notation U V for open subsets of M will mean that U is relatively compact 
in V. 



2.2 The Levy process 

Let us first recall some basic facts about Levy processes At with values in R m . A 
Levy measure is a measure \i on R m \ {0} such that J (|A| 2 A l)p(dX) < 00, and a 
Levy process A f without Brownian part and with Levy measure p, is a process which 
can be written by means of the Levy-Ito representation formula 



A t = Kt+ / XN(ds,dX) + / / XN(ds,dX) (2) 

JO J{\X\<1} JO J{\X\>1} 

for some k G M. m , where the random measure N(dt, dX) = 5(t,AA t ) is a Poisson 
measure on M + x M m with intensity dt[i(dX), and N(dt,dX) = N(dt,dX) — dt fi(dX) 
is the compensated Poisson measure. When J (|A| A l)/i(dX) < 00, then At has finite 
variation, and (j2J) can be simplified as 



A t = Ko t+ [ [ XN(ds,dX). 
Jo ii m 



The relation between k and kq is easily written, and At is a pure jump process 
when k, = 0. We will assume an approximate self-similarity condition and a non- 
degeneracy condition on the Levy process written as follows. 

Assumption 1. There exist some < a < 2 and some positive c and C such that 
cp 2 - a \u\ 2 < [ (X,u) 2 p(dX) < Cp 2 - a \u\ 2 

for u G M m , < p < 1, and where (., .) is the Euclidean inner product. If a = 1, we 
suppose moreover that 



lim sup 

elO 



X p(dX) 

{e<|A|<l} 



< OO. 



If a < 1 (finite variation case), we suppose that A t is pure jump (kq — 0). 
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The process A has finite variation if and only if a < 1. Notice that no smoothness 
is assumed on \i\ it can for instance have a countable support. The additional 
conditions of the cases a = 1 and a < 1 are needed to apply the results of [25, 26J. 

As an example, we can consider stable Levy processes of index a, or more gen- 
erally semi-stable processes (see Chapter 3 of [27J for definitions and properties con- 
cerning these processes). If A t is a semi-stable process and if \x is not supported by 
a strict subspace of IR m , then Assumption [1] is satisfied when 1 < a < 2. When 
< a < 1, the additional conditions mean that A t should be a strictly semi-stable 
process. Actually, Assumption [1] is only concerned with small jumps, so that we can 
also consider for instance truncated semi-stable processes (where jumps greater than 
some value are removed). 

2.3 The equation 

Let us now introduce the process X t , solution of ([T]). When M = R d , the meaning 
of this equation is 

ft rt 



X t = X u + ! a(X s .)dA s + ! b(X s )ds 

J u J u 



\ ( 3 ) 

+ (a(AV,AA s )-X s _-a(X s _)AA s 



u<s<t 



for u < t, where a(x) is the differential at of A i— > a(x,X). Under convenient 
smoothness conditions, the Ito integral with respect to A is well defined, the sum 
converges, and the equation has a unique solution for any initial condition X = xq. 
Notice that the jumps are given by X t = a(X t _, AA t ). Ito's formula enables to write 

f(X t ) = f(x )+ [\Dfa)(X 8 -)dA,+ [\Dfb)(X s )ds 

Jo Jo ^ 

+ J2 ((foa)(X s ^AA s )-f(X s _)-(Dfa)(X s ^)AA^ 

0<s<t 

for smooth functions /, and where Df denotes the differential of /. 

This formula can be used to give a meaning to (JT|) when M is a manifold. We 
consider a coefficient a and a vector field b, 

a : M x R m -»■ M U {oo}, b : M TM. 

We again let a(x) be the differential of A i— >■ a(x, A) at 0, so that a(x) is a linear map 
from R m into T X M: 

a: M -> L(R m ,TM) x ^ d x \ x=Q a(x , X) . (5) 
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We say that (X t ; t > 0) is a solution of (Op) if X t is a cadlag (right continuous with 
left limits) process with values in M U {oo}, which is adapted with respect to the 
completed filtration of At, and such that for any smooth function / and any compact 
subset K of M, Equation (j3J) holds true up to the first exit time from K. We say 
that the process dies at time t if X s G M for s < t and X t = oo. This occurs if 
either Xt- = oo (the solution of the equation explodes at time t), or if X t - G M and 
Xt = a(X t -, AA t ) = oo (the process is killed at a jump of A). 

If we use a proper embedding X : M — > Mr, we can apply (T4j) to the components 
of the embedding X and deduce the equation that should be satisfied on X(M) by 
T(X t ); in order to obtain an equation on Mr , we have to extend the coefficients 
out of Z(M); it is then sufficient to solve the equation in M N and prove that the 
solution remains in M. This is this point of view which is generally used in order to 
prove the existence and uniqueness of a solution, see [H] in the case without killing 
(a(x, A) G M for any x G M). We are here in a slightly different framework (possible 
killing) , and shall prefer to give a proof for the existence and uniqueness by means of 
local charts; this is because the proof based on local charts is then used for studying 
the smoothness of the law of X t (it is clear that the embedded process does not have 
a density in M. N , so proving the existence of a density on M cannot be made through 
a Malliavin calculus on M. N ). The equation on a local chart can also be deduced from 



for / bounded and C°°. In particular, if / is constant on a neighbourhood of x, then 



where \i x is the image of the Levy measure \i by A (->■ a(x, A). This measure describes 
the intensity of jumps; it is the Levy kernel of the process X. 

2.4 Canonical equations 

Up to now, we have been given an equation ([T|) on the manifold M, and have ex- 
plained the rigourous meaning of this equation; in this explanation we need the 



m 



The infinitesimal generator of X t is 




(6) 
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function a given by (jSJ) • It is clear that many functions a are associated to the same 
a. However, there is a particular class of equations, called canonical equations, and 
which were introduced by [21] (see also [IB]), for which a and a are in one-to-one 
correspondence. 

Let us first consider a smooth field a(x) G L(R m , T X M), and let 

a(x,\) = X\(l) for x\(t) = x + / a(x\(s))\ds, (7) 

Jo 

assuming that the solution of this ordinary differential equation does not explode. 
Then it is easily seen that a and a are related to each other by ((SI), and x h-> a(x, A) 
is a diffeomorphism of M onto itself with inverse x i— >■ a(x, — A). Notice that x 
and a(x, A) are in the same connected component, so the study can be reduced to 
connected manifolds. 

For canonical equations, the assumptions needed for our main result (Theorem 
[2] below), or at least a sufficient condition ensuring that they are satisfied, can be 
written in the following simple form. 

Theorem 1. Let A be a Levy process satisfying Assumption^, and let a and b be C°° 
functions on M, with values respectively in L(R m , TM) and TM. Consider Equation 
(pO) with coefficient a given by (jTj), assuming that the ordinary differential equation 
never explodes on M . We suppose that the jumps of A are bounded, and that the 
linear map a(x) : R m — > T X M is surjective for any x; if a < 1 we also suppose that 
6 = 0. Then (jTJ has for any initial condition x a unique solution X t , and the law 
of X t is C°° for any t > 0. 



2.5 Assumptions on the equation 

We now return to the case of a general coefficient a. Let us give the assumptions on 
Equation (CQ) which will imply the existence and uniqueness of a solution, and the 
smoothness of the law of this solution. 

Assumption 2. The conditions on the coefficients a and b are as follows. 

1. Consider, for any e > 0, the map which sends a measure v to the measure 
on M 

{alv){A) = J J l A {a{x,\))l {]M>£} v{dx)fi{d\). (8) 

Let K be any compact subset of M. If e is small enough, then al is a continuous 
map from Ai K into Ai (M) (see Section [Ol for the definition of these spaces). 
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2. Let K be any compact subset of M. There exists e > such that x i— >■ a(x, A) 
is C°° on K for ^-almost any A such that |A| < e. The function a given by ([5]) 
is assumed to exist and to be C°° . Letting D be the differentiation operator on 
M with respect to x, there exists some a V 1 < 7 < 2 such that 



for any C°° real function f , for x G K , for ^-almost any |A| < e, and for any 
J- 

3. The coefficient b is a C°° vector field on M with values in TM. In the case 
a < 1, we suppose that 6 = 0. 

4- The linear map a(x) G L(W m ,T x M) is surjective for any x. 

Here are some comments about these four conditions. 
First condition. This condition states that jumps preserve the smoothness of the 
law of the process. Suppose that x 1— > a(x, A) is a C e diffeomorphism of M onto 
itself with inverse y t— > A). Fix a Riemannian metric on M and the associated 

measure dx\ then the density of a\v is obtained from the density p of v by 

means of the classical formula 



(the determinant is computed for orthonormal bases on the tangent spaces). It is 
therefore sufficient to estimate the derivatives of a -1 ; the condition is for instance 
satisfied for the canonical equations of Theorem [TJ Moreover, the process can be 
killed when it quits some open subset. More precisely, consider a process on a 
manifold Mo associated to an equation with coefficient ao, let M be an open subset 
of M , and kill the process when it quits M. This process is obtained by considering 
the equation with coefficient 



Then a% is continuous if the same property holds true for ao, because the map which 
sends a measure on Mq to its restriction to M is continuous from A4 £ (Mq) to Ai e (M). 
We shall however notice that Assumption [3] below often fails for this example. Other 
examples will be given in Section 12.71 






otherwise. 




(9) 
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Second condition. In this condition the smoothness of a(x, A) is assumed with respect 
to x for A small, but no smoothness is assumed with respect to A except for A — > 0. If 
however a is smooth in (x, A), then the condition is satisfied for 7 = 2. In particular, 
canonical equations of Section 12.41 satisfy the condition. 

Third condition. The additional assumption on b in the finite variation case a < 1 
means that the solution X of our equation is a pure jump process; it is required in 
order to apply the results of [251 EE] . 

Fourth condition. The surjectivity of a(x) is a non degeneracy condition. This 
condition says that small jumps go in all the directions, and is similar to the ellipticity 
condition for continuous diffusions. 

Let us now give the additional assumption concerning big jumps. It is stating, 
roughly speaking, that the process cannot come from a very distant point by jumping. 

Assumption 3. If U is relatively compact, then 



is also relatively compact. 

This assumption is trivially satisfied when M is compact. It is also satisfied in 
Theorem [1] because the jumps of A are in some ball B of K m , so 



and the relative compactness of this set follows from the relative compactness of 
U x B and the continuity of a. The assumption often fails when the process is 
obtained by killing as explained in fl9]); difficulties generally arise when the original 
process can jump from M \ M into M; a counterexample, showing that the law of 
X t is not always smooth in this case, will be given in Section 15.31 

2.6 The results 

The main smoothness result for the solution of ([1]) is the following one. 

Theorem 2. Under Assumptions^ and{^ the equation ([1]) has a unique solution 
X t for any initial condition X = xq, and the law of X t is absolutely continuous for 
any t > 0. If moreover Assumption^ holds true, the law of X t is C°° for t > 0. 
More precisely, if K is a compact subset of M, if t > 0, and if p(t, the 
density of X t with respect to some C°° reference measure with positive density, then 
the derivatives of p with respect to x satisfy 





a~\U) C a(U x B), 



D j p(t,x 0} x)\ < C^ K t-^l' 



(10) 
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uniformly for < t < t , x in M and x in K. 

In the case of canonical equations, we see that Theorem [1] is a corollary of this 
result. In this result jumps of A were supposed to be bounded. When M is compact 
but jumps of A are unbounded, we see that the condition which may cause problems 
is the first condition of Assumption |2j More precisely, if the vector field a(.)A has for 
some A = Ai a stable equilibrium xi, then jumps AA = cAi for large c concentrate 
the mass near xi, so that the density of a%v may be unbounded near x\. 

When the big jumps condition (Assumption E]) is not satisfied, the conclusion 
of the theorem still holds if these big jumps are smooth; some other cases will be 
studied in Section |U and on the other hand, counterexamples will be given in Section 

El 

Theorem 3. Under Assumptions^ and\^ suppose that there exists a decomposition 
/z = // + /z" of the Levy measure such that only ^ satisfies Assumption^ and /i" is 
finite. Let /i| be the image of /J by the map A n- a(x, A); suppose that /J x is C°° , and 
that x i-)- |a4|^ i s bounded on M for any compact K and any I (see the definition 
of this family of semi-norms in Section lKT]) . Then the solution X t of ([T]) has a C°° 
law for t > 0, satisfying (JTUJ). 

The kernel /x| is the part of the Levy kernel \i x coming from fv. The assumption 
of Theorem El therefore requires the part of the Levy kernel for jumps coming from 
distant locations to be smooth, whereas the Levy kernel could be purely atomic in 
Theorem [21 

2.7 About the jump coefficient 

In this article, we study Markov processes X t with infinitesimal generator of the form 
(jSJ), with coefficients a and b satisfying some smoothness assumptions. This covers a 
class of Markov processes, but not all of them. In particular, in this model, the set 
of (times of) jumps of X is contained in the set of jumps of the driving Levy process 
A; the inclusion may be strict, since one may have a(x, A) = x for some A ^ 0, 
but introducing such a behaviour generally destroys the smoothness of a. Thus X 
and A have more or less the same times of jumps, and the rate of jump of X is not 
allowed to depend on its present state x. This is a drawback of this approach, as well 
as of other approaches based on different versions of Malliavin's calculus. A weak 
dependence of the rate of jumps with respect to x can however be obtained through 
Girsanov transformation, see arid cases of more general dependence have been 
studied (under frameworks which are different from ours) by [TH El [15]. We now 
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verify that such a dependence is possible if we drop the assumption of smoothness 
of a, and that this is compatible with our assumptions if only finitely many jumps 
are concerned by this behaviour. 

Suppose that the generator of X is C + where C satisfies the assumptions of 
Theorem [2] or [31 and 



where x h-> h' x {M U {oo}) is finite and C°°. We also have to assume that the kernel 
p! x is Borel measurable, that v i— > v' with 



is continuous from hA l K into Ai (M) (as in Assumption [2]), and that jumps corre- 
sponding to this kernel satisfy Assumption [3] (the extension with the assumptions of 
Theorem [3] is also possible). 

In order to prove that this case enters our framework, we apply the fact that the 
measurable space M U {oo} can be viewed as a Borel subset of K+ (it is a Lusin 
space, see [9]); thus we can view fi' x as a measure on R + , and if we define 



the image of the Lebesgue measure on R + by a(x, .) is the measure n' x on R + , 
plus an infinite mass at +oo; notice in particular that a(x, u) = +oo if u is large 
enough. Let us then introduce a real symmetric Levy process A' t with Levy measure 
n'(d\') = a\\ ! \^ a ^ l d\' , independent of A; then (A, A') is a Levy process on IR m+1 
with Levy measure fi(d\)5 (d\') + 5o(d\)fi'(d\') satisfying Assumption [H We then 
consider the equation driven by (A, A'), with coefficient a(x, A, A'), where a(x, A,0) 
is the coefficient associated to the part £ of the generator, and 



In particular, a(x, 0, A') = x if A' < fi' x (M U {oo})- 1 / The image of // by a(x, 0, .) 
is n' x plus an infinite mass at x, so the solution of the equation has generator C + C 
as required. On the other hand, if K is a compact subset of M and if 





(11) 



a(x, u) = Ml v>0; fi' x ({0,v}) >u \ eR+U {+oo}, 




if A' > and this quantity is finite, 
otherwise. 
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then, with the notations (jSJ) and (TTTT) . 

(a(.,0, .)»(<&) = u'(dx) + - fj,' x (M U {oo}))u(dx). 

Thus these jumps satisfy the first part of Assumption [21 though x i-> a(x, 0, A') is 
generally not continuous. 

We have already seen in (Q that we can consider hard killing of a process (we 
kill it when it hits an obstacle), but this may cause difficulties with Assumption [3j 
With the construction we have just described, we can also consider soft killing where 
the process is killed at some rate h(x) > depending smoothly on x; this means 
that we add the term Cf(x) = —h(x)f(x) to Cf(x), and fjf x is the mass h(x) at oo; 
in this case, the measure v' of (ITT)) is the zero measure on M, so v \- > v' is trivially 
continuous. 



3 Proof of theorems 2 and [3 



The two theorems are proved in several steps. 
3.1 Construction of the solution 

In order to prove the existence and uniqueness of a solution of (JTJ), we first write the 
equation on a local chart (U, $), where $ is a diffeomorphism from an open subset 
U of M onto an open subset V of lR d . We can restrict ourselves to atlases such that 
<3> is the restriction to U of a smooth map on M. If r is the exit time of X from U, 
then (j3J) applied to the components of $ shows that Y t = $(X t ), t < r, should be a 
solution of an equation Q on V with new coefficients 

a*(y, A) = A)), 6*(y) = D$($- l (y))b(<l>- l (y)), 

a*(y)=D$($-\y))a($- 1 (y)). 

More precisely, X is solution of if it is a cadlag process on ¥ U {00} with 
initial condition X = x , satisfying the conditions: 

• For any local chart (U, $) and for any time u, if r is the exit time after u of 
X from U, the process Y t = $(A t ) satisfies the equation ()3j) with coefficients 
(a$, 6$) on {n < t < r}. 

• The jumps of X are given by X t = a(X t -, AA t ). 

• If X t - or Xt is at 00, then X s = 00 for any s > t (00 is a cemetery point). 
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In order to solve the equation (CQ) up to the first exit time from U, we shall have to 
extend coefficients (a$, 6$) out of V, and solve the resulting equation flH]) on M. d . 

Lemma 1. The equation ([I]) has a unique solution X t for any initial condition Xq. 

Proof. Consider open subsets of M with relatively compact inclusions U\ <s U2 <s U3. 
We suppose that there exists a diffeomorphism $ from U3 onto an open subset V3 
of R d , so that (E/ 3 ,$) is a local chart. We define Vi = V" 2 = $(C/ 2 ), and let 

/i : M. d — > [0, 1] be a smooth function such that h = 1 on Vi and /i = on V^. The 
coefficient &$(?/) of (TT2i) is defined on V3. On the other hand, from Assumption [2], 
there exists e > such that 

|A| < e => a{U 2 ,X) C t/ 3 - 

Then a$(y, A) is well defined for y 6 V2 and |A| < e, and takes its values in V3. Thus, 
for y G M d and |A| < e, we can define 



{a${y, A), 6*(j/)) = 




which is an interpolation between (a$,&$) on Vi and the motionless process on the 
complement of V2. Notice that from Assumption [2, 

a*(y,\) =3/ + %)a*(j/)A + 0(|An 

for a$ given by ( IT2|) . and similarly for its derivatives. We can consider on M d the 
equation with coefficients (a$,6$), driven by 

Af = A t -^AA S 1 { |AA S |> £} . (14) 

Our smoothness assumptions on a and b imply that it has a unique solution Y t for 
Y = $(x ) fixed. Let 

r = infjt > 0; Y t £ V x or |AA t | > e\. 

Defining X t = Q~ l (Y t ) for t < r, and X r = a(X T _, AA r ), the process X is solution 
of our equation (p]) up to the time r; if |AA T | > e and X r E Ui, then we can solve 
again the equation from this time r and point X T ; since jumps greater than e are in 
finite number on any finite time interval, we deduce the existence of a solution X t 
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up to its first exit time from U\. Conversely, we can go from X to Y and deduce the 
uniqueness of X from the uniqueness of Y. Thus Equation (JTJ has a unique solution 
up to the first exit time from U\. 

Let us deduce the existence of a solution (X t \ < t < t ) for a fixed t > 0. 
We consider a locally finite atlas of M made of subsets U (k) (any compact subset 
intersects finitely many of these Uo{k)) such that Uo(k) is relatively compact in an 
open subset U(k) of the type of the set U\ of the first part of the proof. We can 
choose for any x an index k(x) such that k is measurable and x G Uo(k(x)), and 
we solve the equation starting from x up to the first exit time T\ of the set U(k(x)) 
(apply the first part of the proof). Let P x be the law of this solution. If 5 is some 
Riemannian distance on M, one can check that sup i<u S(x, X t ) under F x converges 
in probability to as u \. 0, uniformly for x such that k(x) = k (this follows from 
the similar property satisfied by Y = $(X)); thus there exists u k > such that 

P* [n > u k ] > P, [sup 8(x, X t ) < 5(U (k), U(k) c )] > 1/2 (15) 

if k = k(x). 

The equation is then solved by means of the following iterative procedure. For a 
fixed initial condition xq, we solve the equation from time r = up to the exit time 
n from U(k(xo)). If t\ > to (the time up to which we want to solve the equation), we 
have obtained the solution up to time t and we can stop the procedure; if X T1 = oo 
the process stays at oo and the procedure can also be stopped; otherwise, starting 
at time T\ from X n , we solve the equation up to the exit time T2 from U(k(X n )), 
and so on. We stop the procedure either when Tj > to, or when the process has been 
killed (jump to oo). Thus the procedure goes on forever when Tj < to and X Tj ^ oo 
for any j. On the other hand, denoting by Tt the filtration of A t , we have from ( TT5|) 
and the strong Markov property that 

Pfo+i - T d > u k | TV,] > 1/2 on A) = {Tj < to} D {k{X T .) = k}. 

We deduce that 

£ P[AJ] <2^P[^n {r j+1 - Tj > u k }] < 2(1 + to/u k ) 

j>0 j>0 

because there are at most t /u k disjoint intervals of length > u k included in [0,io]- 
Thus, for k fixed, Aj cannot be satisfied infinitely many times. We deduce that if the 
procedure goes on forever, then k(X T ) tends to infinity, so X r _ = oo for r = limrj 
(the solution explodes at time r). In this case the solution is obtained by putting 
X t = oo for t > t. 
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The uniqueness can be proved by considering the first time r at which two solu- 
tions X} and Xt differ, and by supposing that r < oo with positive probability; then 
= X%_, so X\ = X% and the uniqueness of the solution in U(k) for k = k(X^) 
leads to a contradiction. □ 



3.2 The case with only small jumps 

In Lemma Q] we have worked out a construction of the process by means of local 
charts. We now verify that this construction also provides a smoothness result for 
the law on these local charts. In all the proofs, we choose a Riemannian metric on 
M, and the associated C°° measure dx with respect to which we will consider the 
densities of M-valued random variables. We shall study the solution X e of Equation 
(TTJ driven by the Levy process without its jumps greater than e (the process A e 
defined in (fH]l ). The parameter e will be fixed and will be assumed to be small 
enough; notice however that the constants involved in the calculations will not be 
uniform in e; this does not cause any difficulty because we shall never take the limit 
as e 1 0. 

Lemma 2. Consider open subsets of M with relatively compact inclusions Uq (e 
Ux d= Ui (e Us, such that U3 is diffeomorphic to an open subset of M. d . Moreover 
let B be a neighbourhood of the diagonal in M x M. We consider the process X e 
solution of Equation (JTJ driven by A e . The surjectivity ofa(x) is only assumed on 
Ui- Let £ > 0. Then the following properties hold true if e is small enough. 

1. The law of the process X s starting from xq G U\ and killed at the exit from U\ 
has a C l density on Uq. 

2. The density x 1— > q e (t, Xq, x) of this killed process and more generally its deriva- 
tives of order j < £ are uniformly dominated byt~^ dJr '^l a forO < t < t , x G U\ 
and x G Uq. 

3. The density and its derivatives up to order £ are uniformly bounded for < 
t <t and (x , x) G (Ui x U ) \ B. 

Proof. We apply the construction given in the proof of Lemma (TJ denote Vi = ^(U), 
and obtain the process Y t solution of ([3]) driven by A £ , with coefficients (a$, 6$) given 
by f fT3|) . Then X e can be written as X\ = $~ 1 (Y" t ) strictly before the first exit from 
U\, so the killed processes X s and $ -1 (Y) coincide. The non killed process Y t has of 
course not a smooth law for any initial condition, since it is motionless out of Vi- It 
is however possible to modify Y by adding extra independent noise in its equation, 
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without modifying the killed process. For A G R m , |A| < e and A' G M , we can 
replace a$ by 



a$(y, A, A') 

so that 



%)a*(y, A) + (1 - h(y))(y + A') if y G V 2 
y + X' otherwise, 



My, A, A') = y + h{y)a*{y)X + (1 - %))A' + 0{\\[<). 

Then, letting A' be a (^-dimensional Levy process independent of A, satisfying As- 
sumption [1] and with jumps bounded by e, we can solve the equation (jBJ) with coef- 
ficients (a$, 6$) and driven by (A £ , A'). The advantage is that now the differential of 
a$ with respect to (A, A') at (0, 0) is now surjective onto M d , uniformly in y; moreover, 
if e has been chosen small enough and if |A| and |A'| are < e, the map y a$(y, A, A') 
is a diffeomorphism from M. d onto itself, and its Jacobian determinant is uniformly 
positive. Consequently, we can apply Theorem 1 of [25], and deduce that Y t has a 
smooth density y h> py(t,yo,y), with derivatives satisfying 

su P \Dip Y (t,y ,y)\ < C 3 r^' a . (16) 

2/0,2/ 

Moreover, it is proved in Lemma 2 of [26J that py and its derivatives up to order I 
are actually bounded as 1 1 if the number of jumps necessary to go from y to y is 
large enough; thus, for c > fixed and if e has been chosen small enough, 

supj|/>Vy(t,2/o,2/)|; < t < t , \y - y \>cj <Cj. (17) 
If r is the first exit time of Y from Vi, we have from 

E yo [f(Y t )l {t<T} ]=E yo [f(Y t )] -E yo [f(Y t )l {t > T} ] 
and the strong Markov property that the process Y killed at r has density 

Qr(t, 2/o, y) = Pr(t, 2/o, y) - E yo [p Y (t - r, Y T , y)l{t> T }] ■ 



The first term is estimated from ([16]) and (|17p . and for the second one, we notice that 
Y T ^ Vi, and that the number of jumps necessary to go from V° into Vq is large if e 
is small enough, so y py (t — r, F r , y) and its derivatives up to order £ are bounded 
on Vq. We deduce the smoothness of qy, and, by applying the smoothness of 
the law of the process X s killed as well as the estimates claimed in the lemma. □ 
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Lemma 3. Let U be a relatively compact open subset of M. The surjectivity ofa(x) 
is only assumed on the closure ofU. Consider again the solution X s of ([T]) driven by 
A e . If e is small enough, then X £ has a C e density p £ (t, x®, x) on U for any xq G M; 
the density and more generally its derivatives of order j < i are uniformly dominated 
by for < t < t , Xq in M and x in U. 

Proof. It is sufficient to prove the result for U = U , for open subsets U <s U% <e 
U% (s U3 (s C/4 d M such that f/ 4 is diffeomorphic to an open subset of IR d , and 
a(x) is surjective on [/ 4 ; the subset U of the lemma can indeed be covered by a finite 
number of such sets U . Put r = 0, and 

r' k = infjt > r k - XI i U 3 }, r k+l = infjt > r f k ; X e t G U 2 }. 

We can associate to this sequence of stopping times an expansion for the law of Xf 
on U , 

00 

¥ X0 [X £ t G dx] = ^P, [X t £ G dx, r k < t < r' k ] 

k=0 
00 

= E E 

00 
k=0 

where J-"t is the filtration of A and Q e is the transition kernel of the process X s killed 
at the exit from U3. From Lemma [21 this kernel has for any i a C density q £ on Uq 
if e is small enough, so the law of Xf is absolutely continuous on [7 with density 

00 

p £ {t,x ,x) = J^E. 

fc=0 

= q £ (t,x ,x) + J^E. 

k=l 

We already know from Lemma [2] that the first term (which is if xq ^ C/3) and its 
derivatives are dominated by Moreover, we can choose e small enough so 

that the process cannot jump from into U\, and in this case X £ k ^ U\ for k > 1; 
in particular (X £ , x) remains out of a neighbourhood of the diagonal of M x M for 



F xo [X e t edx,t<r' k \ JF Tk ] 1 {t>r k } 



Q £ {t - r k ,X £ k ,dx)l {t > Tk} 



q £ (t-r k ,X £ Tk ,x)l {t > Tk} 



q £ (t-T k ,X £ k ,x)l { t> Tk } 



(18) 
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x G Uq. Thus q £ (t — r k ,X £ k ,x) and its derivatives are uniformly bounded on Uq for 
k > 1 (third assertion of Lemma [2]). Thus the proof of the lemma is complete from 

oo 

\D j p E (t : x Q ,x)\ < Cj(l U3 (xo)r^ a + J2 F *o [t > r k ]) (19) 

fe=l 

as soon as we prove that the series converges and is bounded. We have similarly to 
(TT51) that P x [to > it] > 1/2 for x G ^2 if w is small enough, so by applying the strong 
Markov property of X, 

F xo [r' k >u\ F Th ] >F X0 [r' k -r k >u\ T Tk ] > 1/2 

for k > 1 on {r^ < 00}. Thus 

00 00 

Px,, [X e t eU 3 }>J2 F *o [n < t < r' k ] = E„ 

k=l k=l 
^ 00 

>-J2 V -o[t>r k ] 

k=l 

for t < u, so the series in ffT9l) is bounded, and we have 

|W(t,x ,a;)| < ^(l^o)*-^ 7 " + 2F Xo [X^ G U 3 ]) 
for t < u. If t > u, we use the Markov property, write 

p £ (t, x , x) = K xo [p £ (u, X £ _ u , x)] 

and deduce 

\DY(t,x ,x)\ < C 3 (F X0 [Xl u G U 3 }r^ a + 2F XQ [X £ t G Us]) 

< r (d+j)/a F X0 [X £ t G U A ] (20) 

by using the fact that F[X £ G U^\X^_ U ] > 1/2 on {Xf_ u G U 3 } if u has been chosen 
small enough. □ 

We have proved the estimate (TSUI) which is more precise than the statement of 
the lemma. This property will be used in Section HI 

The absolute continuity of the law of X t claimed in Theorem [2] follows easily from 
Lemma [3] by conditioning on the last jump of A before time t greater than e. If r 



F X0 [r' k >t I T Tk ]l {t > Tk} 
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is this last jump (we put r = when there is no big jump), then we deduce that X 
has a density given by 

p(t, x , x) = E xo [p £ (t - r, X T , x)] . (21) 

However, this formula is not sufficient to obtain the smoothness and even the local 
boundedness of p, because p £ (t — t,X t ,x) is of order [t — r)~ d / a , at least when X T 
and x are close to each other, and (t — r)~ d / a is not integrable if d > a. 

3.3 The case with big jumps 

In Lemma El we have proved the smoothness of the law when big jumps of A have been 
removed. We now have to take into account the effect of these big jumps. Notice 
that the following lemma completes the proof of Theorem [2] when M is compact 
(choose U = M). 

Lemma 4. Let U be a relatively compact open subset of M. There exists a Markov 
process Y t such that the laws of X and Y killed at the exit from U coincide, and Y t 
has on U a C density py(t,yo,y) , the derivatives of which satisfy 

\D i Mt,yo,y)\<c j t-^ a (22) 

for < t < to, yo in M and y in U. 

Proof. For any e > 0, we have a decomposition £ = (£ — £ e ) + £ e of the infinitesimal 
generator, where £ — £ £ is the generator of the process X s driven by A £ , and 

C e f(x)= [ (f(a(x,X))-f(x))fi(d\). 

J{\X\>e} 

Choose U <e U' <e U" <e M. Fix a Riemannian metric on M and consider the 
Riemannian exponential function; then its inverse exp" 1 y £ T X M is well defined and 
smooth if y is close to x. Thus, if e is small enough, we can consider ao(x,\) = 
exp" 1 a(x, A) for x £ U" and |A| < e. Let ho : M — >■ [0, 1] be a smooth function such 
that ho = 1 on U' and ho = on U" c , and let 

a(x,X) = exp x (ho(x)ao(x, A)), b(x) = h (x)b(x) 

for |A| < e. Let £ e be the infinitesimal generator £ — £ £ where coefficients (a, b) have 
been replaced by (a, b); this corresponds to a process X driven by A e , which is inter- 
polated between X s and the motionless process (this is similar to the construction 
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of Lemma CD but we are here on the manifold instead of IR d ); it satisfies Assumption 
El except the surjectivity condition which does not hold out of U" but holds on the 
closure of U'. On the other hand, let h : M — > [0, 1] be a smooth function such that 
h — 1 on U and h = on U' c , and define 

Tf{x) = [ (f(a(x, X))h(a(x, X))h(x) - f(x))»(d\). 

J{\\\>6} 

This means that we consider the jumps A of A greater than e; if the process is at 
a point x before this jump, we kill it with probability 1 — h(x); if it is not killed, 
it jumps to x\ = a(x, A) and is killed with probability 1 — h(xi). We let Y be the 
process with generator £ £ + £ £ ; this is the process X interlaced with jumps described 
by L . This process enters our framework from Section I2.7[ and X and Y coincide 
when killed at the exit from U. 

From Lemma [31 if e is small enough, the process X has on U' a C e density 
p(t, Xq,x) with respect to the Riemannian measure. On the other hand, there are N t 
jumps AA Tfc greater than e on the time interval [0, t], for a random r = (ti, r 2 , . . .); 
we also let To = and append a last TN t +i = t. Let K be the random index k such 
that Tfc + i — Tk is maximal. Then 

^[Y TK+1 - edy\r; A s , < s < r K ] = p(r K+1 - t k , Y tk , y)dy 

on U', with 

\&ftT K+1 -TK,Y n ,y)\<C j (T K+ i-T K y ld+i)/a . (23) 

On {K = N t }, we obtain the conditional density of Y t . Otherwise, we have to 
apply the jump at tk+i to this distribution; we first kill the process with probability 
1 — h(y) and therefore get a C l law on M supported by U'; from Assumption El this 
law is then transformed by al into a C law on M, which is restricted into a C e law 
supported by U' by the second killing. We therefore obtain 

f[Y TK+1 G dy | r; A s , < s < r K ] = p*(y)dy 

for a conditional density p* which is C , with derivatives dominated as in fl23|) . 

This density is then propagated from tk+i to {r K+2 ) — by means of the semigroup 
of X with generator C 6 ; if e is small enough, then x h-> a(x, A) are diffeomorphisms of 
M onto itself for |A| < e, and the process X can be written as X t = $t(Ao) for a flow 
of diffeomorphisms $t of M onto itself; the technique of [12] for compact manifolds 
can be adapted to our case since the process is motionless out of a compact part of 
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M. We choose a copy of $ which is independent of A, and obtain 

P[y TA . +2 _ G dy | r; A s ,0 < s < r K ] 

= Ek(^-^ +1 (y))|det(D^ +a _^ +l )(y) 



(24) 



where E is the expectation only with respect to $, and the determinant is computed 
relatively to orthonormal bases on the tangent spaces. The differential of $ t is 
solution of 

D$ t +dt{x) = Ua($ t {x), dA. t )D$ t (x) + Db{<f> t {x))D$ t {x)dt, D%(x) = I. 

If M is embedded in R , this can be transformed into an equation on M. and 
it is a standard procedure to prove that sup t<T \D<& t (x)\ has bounded moments, 
uniformly in x: use the technique of [13]. By differentiating this equation, the same 
property holds true for higher order derivatives, so actually sup f x \D^ t (x)\, t < t , 
has bounded moments. The same property can be verified for the derivatives of the 
inverse map Q^ 1 , by looking at the equation of its derivative. Thus, in ([2l|), we 
obtain an estimate on the derivatives of the conditional density of Y TK+2 _ similar to 
f[2"31) . By iterating this procedure on all the subsequent jumps r k , we prove that Y t 
has a conditional density; the unconditioned density py is then obtained by taking 
the expectation and satisfy 



\D j PY {y)\ <CE 



on U', because each jump at and each use of the flow $ between and r^+i 
appends a multiplicative constant in the estimation, and they are at most N t of these 
jumps. The number N t of big jumps is a Poisson variable so has finite exponential 
moments. Moreover, 

t 

T K+ i - T K > 



N t + V 
so 

(TJC+I - T K )' (d+j)/a < {N t + l)(d+j)/ at -(d+j)/a 

and we can conclude. □ 

We now give an estimation of the probability for X t to be in some relatively 
compact open subset when one needs many jumps to come from the initial condition. 

Lemma 5. Consider open subsets of M with relatively compact inclusions U n <s 
V n (s U n+ i, and suppose that X cannot jump from U^ +l into V n . Let n > 1. Then 
^xo[Xt £ U ] is 0{t n ) as 1 10 uniformly for x in V^-\- 
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Proof. A similar result was proved in Lemma 1 of [26] for M d ; the proof of this variant 
is much simpler. Let h k , < k < n — 1, be a smooth function with values in [0, 1], 
which is 1 on U k and on V k c . Let C k = sup \Ch k \ for the generator C of X. Then 
h k (x ) — for x G V^_ 1 C V^ c , and Ch k (x) = for a; G £^ +1 , so 

Paat^t e f/ fc ] < E X0 [h k (X t )] = E xo [ Ch k (X s )ds < C k t sup F X0 [X S G U k+1 ]. 

JO 0<s<t 

Applying this inequality for k = 0, 1, . . . , n — 1 completes the proof. □ 

Proof of Theorem® The absolute continuity has already been proved in ( I2T1) . the 
smoothness has been obtained in Lemma H] in the compact case, so we now have to 
study the smoothness of the density in the non compact case. Let Uq be a relatively 
compact open subset of M which is diffeomorphic to an open subset of M. d , and let 
K C Uq be compact. Under Assumption [3j there exists a sequence (U n , V n ) satisfying 
the conditions of Lemma [5] starting from the given Uq. For e and n which will be 
chosen respectively small and large enough, there also exists a process Y constructed 
in Lemma H] for U = U n+ i, with a C density PY(t,yo,y) on U n+ i. The law of Y 
killed at the exit r from U n+ \ has on U the density 

Qr(t,y ,y) = p Y (t,y ,y) - E yo [p Y {t - r, Y T , y)l{ T<t }] ■ (25) 

Consider the process Y when it is out of U n+ i, its big jumps (coming from £ £ ) are 
included in jumps of X, so by construction of the sequence (U n , V n ), it is not possible 
to jump directly from U° +1 into V n ; the small jumps (coming from C £ ) have been 
modified, but they are small, so if e is small enough, it is again not possible to jump 
directly from U° +1 into V n . On the other hand, Y coincides with X on U n+ \ and 
has therefore the same jumps, so it cannot jump from U£ +1 into V k for k < n. Thus 
we can apply Lemma [5] to Y and deduce that P yo [^t G Uq] is 0(t n ) for y V n -\\ 
we also have the uniform estimates (I22p . and we can deduce as in Lemma 2 of 
that if n has been chosen large enough, then Dip Y (t,y ,y) is uniformly bounded 
for < t < to, y G K and y ^ V^_i (this result was proved on R d but U can 
be viewed as a subset of R d ). These estimates on p Y imply by means of (125]) that 
D j q Y (t,y ,y) is 0(t-^ d+j ^ a ) uniformly for < t < t , y G U n+1 , y G K, and is 
bounded if y G U n+1 \ K-i- 

The killed processes X and Y coincide, and the smoothness of the non killed 
process X is then deduced as in the proof of Lemma [3] by considering successive exits 
from U n+ \ and entrances into V n ; we obtain an expansion similar to ([TBI) . When the 
process enters V n , it is at a point out of V n -i, so we have the uniform boundedness 
of the derivatives of q Y starting from this point, and we can proceed and estimate 
the series as in Lemma [31 □ 
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Proof of Theorem^ The decomposition jj, = yU b +yu" of the Levy measure corresponds 
to a decomposition A t = A^+Af of the Levy process into independent Levy processes, 
where A" is of pure jump type. We can apply Theorem [2] and deduce that the process 
X^ driven by A b has a smooth density p b . Let r be the time of the last jump of A* 
before t, with r = when A" has no jump on [0, t]; this last event has probability 
exp — t^iW™). Similarly to ( 12T1) . the density of X t can be written as 

pit, x , x) = E Xo [p (t - t, X T , x)] 

= p\t, xq, x) exp(-t /i fi (M m )) + E xo [p\t - t, X T , x)l {T>0} ] . 

The first term is smooth. On the event {r > 0}, it follows from the assumption of 
the theorem about //jj. that the conditional law of X T given A b has a smooth density 
Pq on M, so 

[p b (t - t, X T , x)l {r>0} ] = E X[ 

We want to verify that the smoothness of po is preserved by p (t — r, ., .). We know 
from the proof of Theorem [2] that p*(t — t,z, .) is C[ on K if the initial condition z 
is out of a large enough subset, so it is sufficient to consider the case of a smooth 
Po with compact support. It also follows from previous proof that it is sufficient to 
prove the result for the modified process Y, and the propagation of the smoothness 
of the law by the semigroup of Y has been obtained in the proof of Lemma HI □ 

4 Lie groups and homogeneous spaces 

In Theorem [21 we have assumed that the process cannot come from infinity by 
jumping (except in the case of smooth jumps of Theorem [3]). This assumption is 
not needed in the affine case M = M d , but other conditions are required in this 
case, and assumptions made for instance in [21] are not intrinsic if M d is viewed as a 
differentiable manifold (they are not invariant by diffeomorphisms) . We can of course 
apply the affine theorem when M is diffeomorphic to M. d , but again the conditions 
on the coefficients will depend on the diffeomorphism. 

Thus, if jumps are not bounded, we need additional structure on M (as this is 
the case on M. d where the affine structure is used). We consider here the case of 
Levy processes on Lie groups; more generally, we consider the case where M is an 
homogeneous space on which a Lie group G acts, and X t is the projection on M of 
a Levy process on G. 



l{r>o} / V (t - t, z, x)p Q {z)dz 
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Exposition about Levy processes on Lie groups can be found in [19] . We shall 
use the theory of integration on Lie groups or homogeneous spaces (Haar measures, 
invariant and relatively invariant measures), which is explained in [6], see also for 
instance [28J. We recall here the points which are useful for our study 



4.1 Lie groups 

Let M = G be a (^-dimensional Lie group with neutral element e and Lie algebra 
g; as a vector space, q is the tangent space T e G; it can be identified to the space of 
left invariant vector fields on G; the Lie bracket of two elements of g, as well as the 
exponential map exp : g — > G can be constructed from this identification. We can 
choose as a smooth reference measure on G a left Haar measure H^L, or a right Haar 
measure %%; each of them is unique modulo a multiplicative constant; they satisfy 

Hl(gA) = nl(A), HliAg- 1 ) = X G (g)nl(A), 
H%{Ag)=H°{A), U%(gA)=XG(g)n%(A), 

for a group homomorphism xg '■ G — > which is the modulus of G. If we are given 
we can define 1-1% by 

H%{dg) = XG(g)nt(dg). (27) 

The group G is said to be unimodular if xg = 1; this holds true when G is compact 
because any group homomorphism from G into must be equal to 1. Let Ad : G — > 
GL(q) be the adjoint representation of G; this means that Ad g is the differential at 
e of the inner automorphism x i— > gxg~ l . Then 

Xa(g)= |detAd fl |. (28) 

On G there are two differential calculi, a left invariant one and a right invariant one. 
The left and right invariant derivatives are the linear forms 



f(g exp(ew)), D^f(g)u = — 

e=0 dE 



/(exp(ew) g) 

=o 



for smooth functions / : G — > R and tiGj. The invar iance means that D^Lhf = 
LhD^f and D^Rhf = RhD^f, with the notations Lhf(g) = f(hg) and Rhf(g) = 
f(gh). The left and right invariant derivatives are related to each other by 

D^f(g)u = D^f(g)Ad- 1 u. (29) 
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If we choose an inner product on g, we can consider the norms \D<_f(g) | and \D^f(g) | 
of the linear forms, and norms corresponding to different inner products are of course 
equivalent. We can also consider classes C[^_ or C%_^ of functions for which the left 
or right invariant derivatives are bounded (without boundedness we can simply use 
the notation C e since the classes for the left and right calculi coincide). 

Let X t be a left Levy process on G with initial condition Xq = e. This is a strong 
Markov process which is invariant under left multiplication, so that its semigroup 
satisfies PtLh = LhPt', equivalently, the infinitesimal generator should satisfy CLh = 
L h C For < s < t, the variable X~ l X t must be independent of (X u ; < u < s) 
and must have the same law as X t _ s . We consider here the subclass of Levy processes 
without Brownian part. Let V be a relatively compact neighbourhood of e in G which 
is diffeomorphic to a neighbourhood U of in g by means of the Lie exponential 
function. Then a left Levy process without Brownian part is characterised by a drift 
k G g and a Levy measure fix on G\ {e} which integrates smooth bounded functions 
/ such that /(e) = D^f(e) = 0; the infinitesimal generator of X can be written in 
the Hunt form as 

£f{g)=D^f(g)K+ [ (f(gx) - f(g) - D^f(g) exp" 1 x) fi x (dx) 

(30) 

+ I {f(gx)-f{g))» x (dx). 
Jv c 

It is explained in [2] that X can be viewed as the solution of an equation driven by 
a Poisson measure on R + x G; by means of a technique similar to Section \2.7\ it is 
also the solution of an equation of type ([1]) driven by a g-valued Levy process. More 
precisely, let ft be the measure on U \ {0} which is the image of (J>x\ v by exp -1 ; on 
the other hand, there exists a bi-measurable bijection i from V c onto a Borel subset 
U' of U c , so we can let \i be on U c the image of ^x\ yc by Then 

Cf(g) = D^f(g)K+ [ (f(g exp A) - f(g) - D^f(g)X)^(dX) 
Ju 

+ [ {figCHXft-fOiDnidX). 

JU' 

Thus X t can be viewed as the solution of X t = a(X t -,dA t ), where A t is a Levy 
process in the vector space g with Levy measure \i supported by U U U', and 



a{g,X) 



g exp A if A G U, 
gi-\X) ifXeU'. 
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Then a(g) : g — >■ T g (G) is the map which sends u to the value at g of the left invariant 
vector field associated to u, so it is bijective, and it is also not difficult to verify the 
second part of Assumption |2j The non degeneracy condition of Assumption [1] on p 
is immediately transferred to an assumption on px as 

cp 2 ~ a \u\ 2 < j (exp^x^Ynxidx) <Cp 2 - a \u\ 2 (31) 

J { exp -1 x\<p} 

for u G Q and p < p small enough so that exp -1 is well defined; the additional 
condition in the case a = 1 is written as 



lim sup 

s4.o 



/ exp 1 x px(dx) < oo. (32) 

J {e<\ exp^ 1 x\<po} 



These two conditions do not depend on the choice of the inner product on g. If 
a < 1, we also suppose that the process is pure jump, so that fl30|) becomes 

£f(a) = / (/W - f{g))f*x(dx). 

Jg 

For the first part of Assumption [21 we have to study the propagation of the smooth- 
ness of the measure by a right translation; this is easy if we choose a right Haar 
measure as a reference measure. 

Thus, under these conditions, we can apply Theorem |2] and deduce that X t has 
a smooth density if px has compact support. If the support is not compact, a 
possibility is to use Theorem [3j Otherwise, we can use the following result. 

Theorem 4. Let X t be a left Levy process on G with Xq = e, the Levy measure px 
of which satisfies ( l3~Tj) . and the additional condition ( |32l) if a = 1. If a < 1 suppose 
moreover that X t is a pure jump process. Then the law of X t , t > 0, is absolutely 
continuous with respect to the left Haar measure H~_. Let i > 0. If 

/ Xg{9)\^ 9 \ 3 fi X (dg) < oo (33) 
Jv c 

for a relatively compact neighbourhood V of e and for j < i, then the density is in 
C*+_. In particular, if 



j \Ad g \ J p x {dg) < oo 
Jv c 



(34) 



for any j , the density is in 
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Proof. The absolute continuity follows from Theorem |2j Assume now (J33l) for j = 
0. Let V be a relatively compact neighbourhood of e, and let fi x and /J^ be the 
restrictions of fix to V and l /c ; as usually, X can be written on [0, t] as a Levy process 
X b with Levy measure fi x interlaced with N t big jumps at times described by fi x . 
We proceed as in the proof of Lemma H] and let r^+i be the end of the longest 
subinterval of [0,t] without big jump. Then, from Theorem [2], conditionally on (t/c), 
the variable Y = X TK+1 _ has with respect to the left Haar measure a C e density py 
satisfying 

\Dlp Y {y)\ < Cj(t k+1 - r K )'^ a (35) 

on V, and the right hand side is integrable as in Lemma HI But Theorem [2] also 
states that the estimate is uniform with respect to the initial condition, so the same 
estimate holds true for any gY, and one deduces that ( 135]) holds uniformly on G. 
Then, we have X t = YZ where Z = X~ _X t is, conditionally on (r^), independent 
of Y. Conditionally on Z and (r^), the variable X t is therefore absolutely continuous 
with density 

p{x\Z) = X g{Z)py{xZ- 1 ) (36) 

(this follows from (l26l)). Thus, by integrating this formula and applying ( l35l) . we de- 
duce that the density of X t is bounded and continuous as soon as Xg(Z) is integrable. 
On the other hand, 

XG{Z)=XG{{Xl K+1 )- 1 Xt) J] XG (X^_X Tk ) 

k=K+l 

where the different terms are conditionally independent given (r^); the process 
Xg(X^) is a geometric Levy process with bounded jumps, so the first term has 
bounded conditional expectation, and we deduce from (133!) for j = that the con- 
ditional expectation of Xg(Z) is bounded by some exponential of N t , so Xg(Z) is 
integrable and the case £ = is proved. 
We can differentiate (1361) and get 

D^p{x\Z)u = X g{Z)D^py{xZ- 1 )M z {u) 

for ii 6 g. For higher order derivatives we have 

Dlp(x\Z)( Ul , ..., Uj ) = X g{Z)DIpy{xZ- 1 ) (Ad z { Ul ), . . . , Ad z { Uj )). 

We deduce the smoothness of the law of X t and the boundedness of its derivatives if 
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we prove that xc{Z)\Kdz\^ is integrable. To this end, we notice that 

Nt 



k=K+l 

N t 



A V,)-^ fc+1 ) 



N t 

n (mk 1 -^)^ } 



k=K+l 



where the different terms are again conditionally independent given (r^); the inte- 
grability is deduced similarly to the case j = by using f[3"3"j) for the terms of the 
second product. The last claim of the theorem follows from (1281) . □ 



Condition ( 1331) is trivially satisfied for j = when G is unimodular, so in this case 
we obtain the existence of a continuous bounded density without any assumption on 
the big jumps. 

Notice also that flMl) is related to the existence of exponential moments for big 
jumps. Assume that fix is the image of a measure fi on g by the Lie exponential (this 
holds for instance when the exponential is surjective). For A G g, let adA : g — > g 
be the adjoint action given by ad A (w) = [A, it] for the Lie bracket [.,.]. We have 
Ad exp(A) = exp(ad A ), so 

|Ad cxp(A) | < exp(|ad A |) < exp(c|A|), 

and (1341) is satisfied if 

exp(C \\\)fi(d\) < oo 



'{|A|>1} 

for any C . This is however the worst case. If G is nilpotent, then the expansion of 
exp(ad A ) is finite, and exponential moments can be replaced by ordinary moments; 
notice however that the class of stable processes introduced by [16] on simply con- 
nected nilpotent groups does not enter our framework, because (|3T|) is not satisfied. 

Theorem H] can of course be translated into the case of right Haar measure, right 
invariant derivatives, and right Levy processes (invariant by right multiplication). 
Then the conditions ( |33l) and ( 1341) are replaced by conditions on Xc{g)~ l and | Ad" 1 1 . 
On the other hand, left and right Levy processes with the same infinitesimal generator 
at e have the same law at any fixed time t, because the right Levy process Y can 
be deduced from the left process X on [0,t] by the formula Y s = X^} s X t . Thus, in 
order to study the law of X t = Y t , one can choose between left and right calculi. 
Notice however that if for instance we apply the result for the right Levy process, 
we obtain that the density p^. of X t with respect to 'K c \ t is of class C^; from ([2T1) . 
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the relation between the densities with respect to the left and right Haar measures 
is p_> = Xg the density is of class C l , but not necessarily of class C{^_ (the 
relation between left and right invariant derivatives is given in (1291) ). 

4.2 Homogeneous spaces 

We now consider the case where the manifold M is an homogeneous space M = G/H. 
More precisely, G is a m-dimensional Lie group of transformations acting transitively 
and smoothly on the left on M, and H = {h G G; h(o) = o} is the isotropy group of 
some fixed point o of M; the projection tt : G — > M is given by ir(g) = g{o). We can 
choose a Lebesgue measurable section S of 7r _1 (which exists from the measurable 
section theorem), and any g in G can then be uniquely written as g = S x h for 
x = Tr(g) = g(o) and some h in H. The action of G on M can be written as 
g(y) = ir(gSy). We will denote by g and \) the Lie algebras of G and H. 

We can look for a measure on M which would be invariant under the action of 
G, but such a measure does not always exist. We therefore weaken the invariance 
into a relative invariance property, see [6j [28] for an introduction to the topic and 
some of the properties which are given below; we say that a Radon non identically 
zero measure % M is relatively invariant under the action of G with multiplier x if 
X '■ G — > is a group homomorphism and 

u M {g{A)) = x (g)n M (A). 

Then the measure is invariant if x = 1 ; this is necessarily the case when G is compact. 
For instance, we have seen in (|26|) that the right Haar measure on G is relatively 
invariant under the left multiplication with multiplier xg- A relatively invariant 
measure also does not always exist, but it exists in more general situations than 
invariant measures. It exists on M = G/H if and only if Xg/Xh '■ H — > IR+ can 
be extended to a group homomorphism \ : G — > W^; in this case there exist a 
relatively invariant measure with multiplier x, and this measure is unique modulo 
a multiplicative constant. In particular, an invariant measure exists if and only if 
Xg — Xh on H; this property holds in particular when H is compact. 

The relationship between left Haar measures on G and H and relatively invariant 
measures on M is the following one. If is a left Haar measure on H and if 
X '■ G — > is a group homomorphism, then a Radon measure "H M on M is relatively 
invariant with multiplier x if an d only if the measure H^L defined on G by 

/ f{9)Ut{dg) =\\ f -{S x h)U M {dx)Ul{dh). (37) 
Jg J Jmxh X 
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is a left Haar measure. For right Haar measures, (137]) becomes (use ( |2Tj) ) 

/ f{g)U%{dg)= [[ f(S x h)^(S a )n M (dx)H*(dh). (38) 
Jg JJmxh X 

If H is a G-valued variable with densities p<_ and with respect to and 
"H^, we deduce from (1371) and (1381) that 7r(S) has density 

p( x ) = / tz(S x h)Hf.{dh) = ^(S x ) [ p^{S x h)H*{dh) (39) 

with respect to "H A/ . Estimates on p can therefore be deduced from estimates on p^ 
or p^, but this is clearly simpler when H is compact. 

We now explain how left and right invariant differential calculi on G can be 
transported to M. For the left invariant calculus, notice that if F is a smooth 
function on M, we can differentiate / = F o it which is defined on G and consider 

D^F{x)u = D^f{S x )u (40) 

for u G g; notice that the result is if u is in f), so the differential is actually a linear 
form on the vector space $j/f). The problem is that it depends on the choice of the 
section S; if S' is another section and if we fix x, then S' x = S x h for some h G H, 
and 

D^f(S' x )u = D^f(S x )M h u. 

If however H is compact, then we can choose on q an Ad (if) -invariant inner product 
(this means that jAd^wl = \u\ for h in H), we can consider D^F(x) on the orthogonal 
p of f) in g, and the norm \D<_F(x) \ will not depend on the choice of S. It is invariant 
under the action of G in the sense 

\DUFog)(x)\ = \(D^F)(g(x))\. 

Higher order derivatives have a similar behaviour, and we can consider the classes of 
functions C[^_ on M. 

For the right differential calculus, we only consider the case where (M, ■) is itself 
a Lie group, and left translations of M form a normal subgroup of G. In this case, 
a canonical choice for the section S is to let S x be the left translation by x, so that 
S : M — » G is an injective group homomorphism. The group G is a semi-direct 
product 

G = S{M) x H (41) 
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satisfying the commutation property hS x = Sh( x )h for h G H. All elements of G can 
be written in the form S x h for some x G M and h G H, and we have the product 
rule 

S(xi)h 1 S(x 2 )h 2 = S(xi ■ h 1 (x 2 ))h 1 h 2 . 

The vector space q can be written as q = f) © m for the Lie algebras f) and m of H 
and M ~ 5*(M). We can consider the differential D^F(x)u computed on the Lie 
group M for u G m, and therefore the classes of functions C^. The relation with 
the differential on G is given by 

D^F(x)u = D^f(S x )u for / = F o vr 

similar to (jlDjl . because exp^ii)^ = 5 , exp ( £tt ) a ,. The behaviour under the action of G 
is 

D^(F o = ( J D^F)(( ? (x))Ad,n. (42) 

If "H^. and are right Haar measures on H and M, a right Haar measure can be 
defined on G by 

/ f{g)U%{dg) =\\ f(S x h)H™(dx)UZ(dh) (43) 
Jg JJmxh 

because the right hand side is invariant under the right action of H and S(M) (use 
S x hS y = S x .h( y )h). By comparing with (138]) . we see that 'K 1 ^ is a measure on M with 
multiplier 

X (S x h) = XG (S x h)/ X H(h), (44) 

since this is a group homomorphism and Xg/x{S x ) = 1- In particular, the Haar 
modulus of M is Xm{ x ) — x(S x ) = Xg{S x ) and (I3"5j) becomes 

p(s) = / p^(S x h)H^(dh). (45) 

In the next subsections, we study some classes of Markov processes on M. In the 
first case, we assume that H is compact, and consider the class of Markov processes on 
M which are invariant under the action of G; they can be written as the projection of 
left Levy processes on G, and we apply the left differential calculus. We also consider 
the case where X is the projection of a right Levy process. In the second case, we 
assume that G is a semi-direct product of type (T4TI) . and we let X be the projection 
of a right Levy process on G; we then apply the right differential calculus. 
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4.3 Case 1: compact isotropy subgroup 

We here assume that M = G/H for a compact Lie subgroup if of G; in particular, 
there is on if a unique probability measure H H which is both a left and right Haar 
measure, and there is on M a measure 7f M which is invariant under the action of G, 
and which is related to a left Haar measure on G by ( |37l) with \ = 1- We choose an 
Ad(ff)-invariant inner product on g; it induces a if-invariant inner product on the 
tangent space T a M, and a G-invariant Riemannian metric on M. We consider the 
class of Markov processes X t which are invariant under the action of G; this means 
that Pt(F o g) — (PfF) o g for any g G G; in particular, the law of (Xt, t > 0) with 
initial condition X = o is invariant under the action of H. From Theorem 2.2 of 
[19] . these processes are obtained as X t = 7r(H t ), where E t is a left Levy process 
on G which is invariant under the right action of if. The invariance of X implies 
in particular that its Levy kernel fi x can be deduced from \x a which will be simply 
called the Levy measure fix of X; it is a if -invariant measure on M \ {o} which 
integrates 5 2 (o,x) A 1 (for the Riemannian distance 5). The Levy process H on G 
can be obtained by taking a Levy measure given by 

fjLs(A) = [ [ l A (hS x h~ x )H H (dh)y. x (dx), (46) 

J JHxM 

where the if-invariance of fix implies that /z= does not depend on the choice of the 
section S (see [IS]). Then the generator of H can be written by means of fl30|) . for a 
neighbourhood Vq = exp Uq, where Uq is a small enough ball of q (in particular Uq 
is Ad(if)-invariant), and for a drift k G g which is Ad(if)-invariant. 

For < s < t fixed, notice that S t has the same law as H S H£_ S , for S' an 
independent copy of S; moreover, the variable S s can be written as Sx s h s , for /i s a 
if- valued variable. Thus, letting X' = 7r(H'), 

X, = H t (o) ~ S s H;_ s (o) = SxMXU ~ 5x s (A;_ s ) (47) 

because the law of A t '_ s is if-invariant. More generally, if fix and /i2 are two H- 
invariant measures on M, one can define the convolution fi\ * fi2 as the image of 
the product measure by (x, y) h> S x (y); it does not depend on S and is again if- 
invariant. The relation ( I47|) shows that the law z/ t of X t satisfies u t = v s * v t - s . If the 
convolution product is commutative on the set of if-invariant measures, then (G, if) 
is said to be a Gelfand pair. 

An example of space M is the hyperbolic space viewed as a subspace of the 
Minkowski space, namely 

d 

B. d = |s = {x ,xi, . . . ,x d ) G R 1,d ; \x\ = 1, x > o|, \x\ 2 = x 2 Q - ^x 2 , 

i=X 
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with o = (1,0, ...,0). It can be viewed as G/H, where G = SO + (l,d) is the 
restricted Lorentz group of linear transformations of M. 1,d which preserve the pseudo- 
norm, the time direction and the space orientation, and H — {h G G; h(o) = o} ~ 
SO(d). Then it is known that (G,H) is a Gelfand pair. 
For our result, we need the functions 

Xm(x) = Xg(S x ), A(x) = \ Ad Sx \ 

which do not depend on S because xg — 1 on H and the inner product of q is 
Ad(if)-invariant. 

Theorem 5. On M = G/H for H compact, let X t be a G-invariant Markov process 
without Brownian part, with X = o, the Levy measure fix of which satisfies the non 
degeneracy assumption fl3T|) where exp -1 = exp" 1 denotes the inverse Riemannian 
exponential function based at o, and the additional condition fl32|) if a — 1. If a < 1 
suppose moreover that X t is a pure jump process. Then the law of X t , t > 0, is 
absolutely continuous with respect to the G-invariant measure 1-L M . Let I > 0. If 

X M {x)A{x) j fi x {dx) < oo (48) 

for a relatively compact neighbourhood V of e and for j < i, then the density is in 
C[^_ (see Section \4~J\ for the definition). In particular, if 

I A{x) j [ix{dx) < oo (49) 
Jv c 

for any j, the density is in C^_. If(G,H) is a Gelfand pair, conditions fj48|) or fj49|) 
are not needed. 

Proof. We write X as X t = vr(SJ) for a left Levy process Bf with Levy measure 
[131 given by (1151) ; denote by CJ its infinitesimal generator. The question is to know 
whether the non degeneracy condition ( 13T1) for X can be translated into the similar 
condition for H' (and similarly for (132|) when a = 1). Recall that q is written as the 
orthogonal sum of f) and p, and we can choose the section S such that S x G expp 
for x in a neighbourhood of o; then S is uniquely determined and smooth on a 
maybe smaller neighbourhood; the measure fix on a neighbourhood of o is therefore 
transported to a measure [i! on a neighbourhood of in p, and // satisfies Assumption 
[H If S x = exp A, then 

exp-^hSvhr 1 ) = Ad h X + 0(|A| 2 ). 



\ 
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On the other hand, we have from (1461) that 



I(p) = / <exp 1 £,u) 2 fi s ,(d£) 



J{\ex P '^\<p} 




{lexp-l^h" 1 )!^} 



(exp 1 {hS x h 



), u^ 2 'H H (dh)fix{dx) 



so obtaining the lower and upper bounds (I3TT) for I(p) is equivalent to estimating 



The upper bound follows easily since p' satisfies Assumption [TJ and for the lower 
bound, we can restrict the domain of integration H to a small neighbourhood of the 
unity on which |Ad^A — A| < c|A| for c arbitrarily small. We deduce that f l3Tj) holds 
true for the process 5' but the lower bound is only for u in p. In order to obtain all 
of q, we add extra independent noise in Ef. Let E" be a left Levy process on H with 
generator C" satisfying the conditions of Theorem HJ it is therefore associated to a 
measure p" on a neighbourhood of in f), satisfying ( 13TI) ; we have C"(F o ir) = 
for any smooth F on M. Let S be the process with generator C = £' + Then 
C(F o 7r) = o 7t), so = vr(Sj) can also be written as X t = vr(S t ), and S t is a 
left Levy process associated to the measure p' + p" which now satisfies (|3T1) for any 
u in q. Condition ( |32l) is similarly extended to 5 when a = 1. 

Thus we deduce from Theorem @] that E t has a density p^(t, .) with respect to 
We have from (|39|) that the density of X t with respect to T-L M is given by 



Condition ( 148|) for A implies ( l33l) for S' because Xc(h) = \Adh\ = 1 for h in if, so 
that XcihSxh^ 1 ) = Xm(%) an d \Ad^s x h- 1 )\ = and ( 1331) always holds for S" for 

the same reason. Thus Theorem H] can also be applied to S for the smoothness of 
the law. In the case I = 0, the continuity of p follows from the continuity of p<- and 
the fact that we can choose a S which is smooth in a neighbourhood of x. In the 
case £ = 1, we have 



from the left invariance of H H , and we deduce from the definition (HOl) that 
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\D^p{t,x)\< / \D^{t,S x h)\n H (dh). 
Jh 

The study of higher order derivatives is similar. 

If (G, H) is a Gelfand pair, we write a decomposition fix = Hx + f^x where 
fi x is the restriction to a i^-invariant relatively compact neighbourhood of o; this 
corresponds to a decomposition £ = C) + C\ and X can be viewed as the process 
with generator C) interlaced with big jumps described by CK Conditionally on the 
times of the N t big jumps before t, the law of X t is therefore the convolution of 
(2N t + 1) i7-invariant laws. From the commutativity of the convolution, all the big 
jumps can be put together, and we can write X t in law for t fixed as S x i{X\\ The 

law of X\ is in C^_, and this smoothness is preserved under the action of S x %. □ 

If (G, H) is a Gelfand pair, the technique of previous proof for putting together big 
jumps can be extended to other cases. For instance, we can obtain the smoothness of 
the law if the Levy measure is the sum of two measures, and only one of them satisfies 
the assumptions. This implies that the upper bound in (jHTl) can be weakened, as in 

We can also consider the class of processes X t = 7r(H t ), for right Levy processes 
on G. There are Markov processes with semigroup P t f(x) = E/(H t (x)). Noticing 
that Ht can also be viewed at fixed time as the value of a left Levy process, we 
again apply left differential calculus and immediately obtain from the above proof 
the following result. 

Theorem 6. On M = G/H for H compact, let X t = 7r(S$) = S$(o) for a right Levy 
process H t on G. We suppose that the left Levy process having the same generator 
at e as H satisfies the assumptions of Theorem^ for some i. Then X t has a C[<_ 
density with respect to the G-invariant measure T-L M . 

4.4 Case 2: Semi-direct product 

We consider as in Theorem O processes X t = 7r(5 t ) = S t (o) where S 4 is a right Levy 
process on G, but do not assume that H is compact. Instead, we suppose that G is a 
semi-direct product as described in ( HT|) . and we apply the right differential calculus; 
recall that the right Haar measure 7£™ of M is relatively invariant under the action 
of G with multiplier \ given by ( I44jl . 

A typical example is when M is the additive group M. d and G is the afline group. 
Then H = GL(d). 
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Theorem 7. On M = G/H for G = S(M) x H, let X t = 7r(H t ) = E t (o) for a right 
Levy process S t on G. We suppose that the Levy measure /i~ ofE satisfies (I3"T1) . and 
the additional condition (1321) if a = 1, and that E is a pure jump process if a < 1. 
Assume also that 

' ' Ad' 1 \ j fi S (dg) < oo (50) 



v- x(g) 

for j < where V is a relatively compact neighbourhood of e, and \ is given by 
Then X t has a C\_^ density with respect to the right Haar measure H^. 

Proof. Let us consider S t . We know that it can be viewed as the solution of an 
equation driven by a q- valued Levy process A. We denote by H e or H(e) the same 
process when jumps of A greater than some e (for some norm) have been removed, 
and by X s or X(e) its projection on M. We have seen in (1201) that we have an 
estimate for the density of S e and its derivatives at e involving P[Hf G Vq] for a 
relatively compact neighbourhood Vq of e, uniformly in the initial condition. From 
the right invariance of the process, we also have an estimate for the density at g 
involving P[Hf G Vag], uniform in g. In particular, we can write 

\DUv% E) {t,S x hg» l )\ < C 3 t-^l a n~l e VgSMo 1 ]- (51) 

for any go G G, x G M, h G H . We want to estimate the integral of this quantity 
with respect to h. Let Vm be a relatively compact neighbourhood of o; then V G = 
Vq SiVu)" 1 is a relatively compact neighbourhood of e. For x in M, let V M = 
x" 1 ■ Vm ■ x. We have Vq C VgS^S'^S'" 1 for y G V^-, so 

P[S? G Kg^H- 1 ] < P[~f G V G S x S y hg Q l ] 

for any y G V M and h & H. By taking the mean value on V^, we obtain 

P[S« G IfeS^" 1 ] < uM<v*\ f PK G WW*), 

and -H^(V^) = "H^ (Va/)/x a/ (x). Thus 

CxaKx) / / P[~f G V^Syhg^UZ {dy)U*{dh) 
C X m(x) J P[H| G V^gg^Ti^dg) 

^XM(x)E[^(5- 1 (yG)- 1 ^o)] = CU^iV')- 1 ) (52) 



< 
< 
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where we have used (1431) in the second inequality and xg^S' 1 ) = 1/xg(S x ) = 
1/Xm(x) hi the last equality. The density of Hfg at g is p^ e '(t, gg$ ), and simi- 
larly for its right invariant derivatives. The law of the variable 7r(2fgo) is the law of 
Xf with initial condition xq = ir(g ), and fH5|) becomes 

p^ e \t,x ,x)= f p^(t,S x hg^)n^(dh). 
Jh 

By differentiating, we have 

\D^\t,x ,x)\ < [ \Di,p% s \t,S x h 9 o 1 )\H^{dh) KCjt-V+M* 
Jh 

from (151I) and ( l52l) . Thus the smoothness of the law of X s is proved. We now 
have to take into account big jumps with the technique of Lemma HJ by considering 
X e interlaced with the big jumps, and letting (t k ,t k+ i) be the longest interval 
without small jumps. Previous argument shows that, conditionally on the times of 
big jumps, the variable X TK+1 _ has a density which is in C$_+, and its derivatives 
are of order (tk+i — r^) - ^- 7 )/". The variable X t is then obtained from the action 
of T = H t H^ , so its density is 

^(t,x ,x)=E[^(T- 1 (x))/ X (T)]. 

The variable l/x(T) is conditionally integrable (given the times of big jumps) if fj50|) 
holds for j = 0, so the theorem can be proved for i = by the technique of Lemma 
HI the case i > 1 is similar by applying ( ]42i) for the derivatives of p^(T _1 (x)). □ 

5 Examples 

We here give some examples, and also some counterexamples where the "big jumps" 
condition (Assumption [3]) does not hold, and the smoothness of the density fails. 

5.1 Isotropic jumps 

We have assumed that X is solution of an equation driven by some Levy process 
A, but usually, jumps are often described by the Levy kernel u, x , the image of \x by 
A i — y a(x, A). It is not easy to know when some Markov process with some Levy kernel 
can be represented as the solution of an equation of our type. We have already seen in 
Section l2TTl how it is possible to deal with a finite part of the Levy kernel. We now give 
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the example of p] where this is globally possible. Let M be a complete Riemannian 
manifold, and suppose that fi x is the image by (r,u) >->■ exp a ,(rw) of £g> v x , where 
fin is a measure on (0, oo) (radial part), and u x is the uniform probability measure 
on the unit sphere of T X M (angular part); this means that we choose a direction 
uniformly in the unit sphere, then go along a geodesic in that direction, at a distance 
chosen according to jiR. Such a \i x is singular if /ir is singular. 

In order to construct an equation for this process, as explained in pQ, we lift it 
to the bundle 0(M) of orthonormal frames, as this is classically done in the Eells- 
Elworthy-Malliavin construction of the Brownian motion. Points of this bundle can 
be written as £ = (x, g) for x G M and g : M d — > T X M is an orthogonal linear map; 
we put 7r(£) = x. Then, for A G M. d , we can define a(£, A) for £ = (x, g) by 

A)) = exp x (g\), 

and the frame at 7r(a(£, A)) is deduced from the frame g at x by parallel translation 
along the geodesic (exp x (g\t); < t < l). Let 5 be the solution of the equation on 
0(M) with this coefficient a and with 6 = 0, driven by a symmetric Levy process A 
with Levy measure /i = \lr <g) z/, where z/ is the uniform measure on the unit sphere 
of M. d . Then X = 7r(S) is the process that we are looking for. The process 5 can be 
viewed as the horizontal process above X. Notice that if (e;) is the canonical basis 
of M. d , the vector fields a(.)ej are the canonical horizontal vector fields on O(M), and 
the equation of H is a canonical equation, since a is obtained from a by means of (13). 

However, the surjectivity of a(£) cannot be satisfied, since the dimension of O(M) 
which is d(d + l)/2 is greater than the dimension d of M d . Nevertheless, we can 
add extra noise in H without modifying the law of X and get this non degeneracy 
condition; more precisely, the extra noise acts vertically on the process 5. We enlarge 
the space M d of the Levy process as M d x 0(d), put a(x, A, e) = a(x,X), and let 
a(x,0,g) be the vertical transformation which modifies the frame by composing it 
with g. We let A' be an independent Levy process on 0(d), and consider the equation 
driven by (A, A'). 

We choose the reference measure on O(M) which, when projected on M is the 
Riemannian measure, and which is on each fibre the uniform measure (normalised 
measure invariant under the action of 0(d)). We consider the process 5 on 0(M) 
the initial condition So of which has uniform law on the fibre above some xq. We 
deduce from Theorem [2] the smoothness of the law of H t on 0(M) if [Ir has bounded 
support, and L\\\< p \ |A|V.r(^) i s bounded below and above by constants times p 2 ~ Q 
as p J, 0. Conditionally on X t = vr(S t ), the variable Et has uniform law on the fibre 
above X t , so the density of H t is a function of x only; this is also the density of X t , 
so X t has a smooth law too. 
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When M is the sphere or the hyperbolic space, then M is a Riemannian symmetric 
space, and we are in the framework of a G-invariant process on M = G/H, where 
(G, H) is a Gelfand pair. Thus the smoothness of the law holds on M. d also in the 
case of unbounded jumps (Theorem [5]). The hyperbolic space is diffeomorphic to 
M. d ; however, if we use normal coordinates, we cannot apply the theorem of [23] for 
M, d because the ellipticity condition is not satisfied. The matrix (oo*)" 1 , which is 
the hyperbolic Riemannian metric, explodes indeed exponentially fast at infinity, 
whereas at most polynomial growth was assumed in [24J. 

5.2 Affine transformations 

Let G be the group of affine transformations of M. d ; it enters the framework of ( 14TT) 
as G = M. d x GL(d), where g = (g 2 , 9i) is the map g(x) = g\x + g 2 . The Lie algebra 
is the set M. d x gl(d), and if u — (u 2 ,ui), we have 

gexp(eu)g~\x) = x + e(g 1 u 1 g~ 1 x + g^u 2 - giUxg~ l g 2 )j + 0(e 2 ), 

so 

Xc(g) = |det#i|, AdgU = (gi^g^ 1 , g^ - giu 1 g^ 1 g 2 )) 

We deduce from Theorem H] the uniform smoothness of the laws of a left Levy process 
H( and its projection X t = 7r(H t ) on Mr, under the non degeneracy condition on the 
small jumps (with additional conditions if a < 1), and if the moments of |Ad s | for 
the big jumps part are finite. Under the similar condition on | Ad" 1 1 , we can also 
consider the right Levy process, and deduce the smoothness of its M. d component 
from Theorem [3 in this case H = GL(d) is unimodular, so x = Xg- 

Notice that G is not unimodular, so without the assumption on big jumps, we 
are not even sure of the local boundedness of the density of X t , except when we can 
apply Theorem [3] (smooth Levy measure). In order to find a counterexample, we are 
going, for simplicity, to consider a subgroup of G. 

So let now G be the group of transformations of IR d generated by translations 
and by dilations of rate e n , neZ. Thus G = IR d x Z, where (y,n) corresponds to 
the transformation z 1— > e n z + y. The composition law of this group is 

(y2,n 2 ).(yi,n 1 ) = {e n2 y x +y 2 ,n 1 +n 2 ). 

Its Lie algebra is Abelian and can be identified to M. d . Denoting by dn and dy the 
counting measure on Z and the Lebesgue measure on R d , the Haar measures on G, 
the modulus and the adjoint representation are 

T-L^_(dy dn) = e~ nd dydn, %z>(dydn) = dy dn, XG(y, n ) — e " d > ^{ y ,n)U = e n u. 
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We consider on G left Levy processes X t = (Y t , N t ) such that N t is a random 
walk on Z, and Y t is deduced from a Levy process At on M d , independent of N, by 
means of ^ 

K= / e N 'dA„. 



We want to study the density near (0, 0) of X\ with respect to or equivalently 
the density near of Y\ restricted to {N\ = 0}. Suppose that A is an isotropic stable 
process. Then 

Yi~ a( [ e aNs ds) ~ ( f e a7Vs dsV /a Ai, 



and its conditional density at given A" is 

p Y (0\N) = c A (l e aN °ds) ' (53) 



where ca is the density of Ai at 0. Let p n be the mass of the Levy measure of A" 
at n, with YlPn < oo. Let A n , n > 1, be the following event: the process A^ has 
exactly two jumps before time 1, one jump of size — n followed by a jump of size +n. 
Then the probability of A n is CNP-nPn with cn = \ exp — ^2pj] conditionally on A n , 
the times of the two jumps are obtained from two independent variables uniformly 
distributed on [0, 1], and by letting the negative jump be the smallest one, and the 
positive jump the largest one. The law of X\ restricted to A n is smooth; by applying 
(j551) and by denoting c = c\c^, its density at (0, 0) is 



p An (0,0) = 2cp^ nPn J J (s + (t - s)e~ na + 1 - tj ^dsdt 
= 2cp^ nPn / (1 - s) 1 + s(e- na - 1) ds 



> cp_ nPn e- 2na (l - (1 - e""") 2 ) ^ 

> cp^ nPn e n{d ~ 2a) . 

The law of X± restricted to Ui<fc<n^fc is smooth with density Y^k=iVA k , so for any 
neighbourhood U of (0, 0), the density p of X\ satisfies 

n 

ess sup p > p_ k p k e k( - d ~ 2a \ 
u t^i 

If the series diverges, then X\ does not have a locally bounded density, and examples 
can be constructed as soon as d > 2a. 
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Choose for instance 

p n = e~ Pn , p. n = e-™ (54) 

for n > 1, /3 > 0, a > 0. In this case the density is not locally bounded if d > 
2a + (3 + a. On the other hand, Theorem H] implies that if d < /3, the density is 
bounded uniformly with respect to the initial condition. 

Consider now the law of Y\. We can consider the right Levy process X' with the 
same generator at e as X; if N and A are the same independent random walk and 
stable isotropic processes, X' can be obtained as X' = (Y', N), where Y' the solution 
of Y t ' +dt = a(Yt, dN t , dA t ), and a(y,n,X) = e n y + A. We have the equalities in law 
X[ ~ X x and Y{ ~ Y x . 

The study of the law of Y' has been the subject of Theorem [71 For the example 
( 1541) . this theorem implies that Y[ ~ Y\ has a bounded density if d < a. On the other 
hand, we have just verified that this density is not locally bounded if d > 2a + (3 + cr; 
actually, this condition can be improved because we now study Yi instead of X\. 
The event A n can be replaced by A' n : the process N has exactly one jump of size — n 
before time 1, and we deduce that the density is not locally bounded if d > a + a. 



5.3 Killed processes 

Consider a process which satisfies our sufficient conditions for the smoothness of the 
law in some manifold, say the affine space R d , but which is killed at the exit from 
some open subset M. The existence of a locally bounded density is preserved, but it 
appears that this killing can destroy the smoothness of this density. We have seen in 
(GO) that Assumption [2] is preserved, but Assumption [3] may fail; a problem can occur 
when jumps are allowed to enter the subset (coming from the infinity of M), or at 
least (by applying Theorem [3]) if the non smooth part of these jumps are allowed to 
enter the subset. If directions of non smooth jumps lie in some closed cone of the 
vector space R d , then smoothness is preserved if the obstacle M c is such a cone based 
at some point of lR d ; in the general case however, roughly speaking, the obstacle can 
produce singularity behind it. Let us give a one-dimensional example. 

Let A t = A° — Aj where A° and A^ are real, independent, and are respectively a 
symmetric stable Levy process with Levy measure |A| _a_1 ciA and a standard Poisson 
process. Let X t be the process At killed when it quits M = (— oo, 1). It is solution 
of the equation corresponding to b = and a(x, A) = x + A if x + A < 1, equal to oo 
otherwise; the initial condition is xq = 0. 

Notice that the similar process without Aj has smooth densities from Theorem 
[31 If A] is added instead of subtracted, the law should also be smooth since the 
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appended jumps do not enable the process to enter M. We are now going to prove 
that in our framework, the law of X\ is not C 1 as soon as a < 2/3. 

Let t be the lifetime of X, which is the first exit time of A from M. Then the 
density of Xi is 

p(x) = q(l,x)-E[q(l-r,x-A T )l {T<1} ] (55) 
where q(t, .) is the density of A t . We have 

q(t, x) = e~* ^2 IT x + k ) 

k 

where the density go of A is C°° on (R + xl)\{(0, 0)}, with bounded derivatives out 
of any neighbourhood of (0, 0). Consequently, q is smooth on (R + xK) \ ({0} x Z_), 
and 

_ t k 

q {k) {t, x) = q{t, x) - e -nQo(t, x + k) (56) 

is smooth on M + x (—k — 1, —k + 1). We are going to study the derivative of p(x) as 
x | 0. 

We need some information on the joint law of the exit time r and the overshoot 
A T — 1. To this end, we first check that r is almost surely a time of jump of A; this 
means that the process cannot creep upward, see [10J; a simple way of verifying this 
fact is to notice that A cannot creep both upward and downward because it has 
no Brownian part; since it is moreover symmetric, it creeps neither upward, neither 
downward, and A which is obtained by adding a process with finitely many jumps 
satisfies the same property. Consider the joint law of (r, A T _, A T ), and denote by T 
the set of jumps of A; notice that r is a jump of A ; then 

E [f(r, A T _, A r )l {T<oo} ] = E /(*> Xt-,X t . + AA f )l { x t _^oo}l{AA t >i-x t _} 

teT 

poo poo 

= E / / f^X^Xt + X^x^yl^.x^X-^dXdt 
Jo Jo 

from a key formula of stochastic calculus on Poisson spaces. Thus 

P[r e dt, A T _ G dx, A T G dy] = P[X t e dx] (y - xY^dtdy 
on (0, oo) x (-co, 1) x (1, oo), and the density of (r, A r ) is 

C(t,y) = E[(y-X t )- 1 - a ] (57) 
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on (0,oo) x (l,oo); in particular y \-t ((t,y) is C£° on [1 + e,oo) for any e > 0, 
uniformly in t. Let h : R — >■ [0, 1] be a smooth function which is 1 on (—00, 4/3] and 
on [5/3, +00). From (1331) we have 

/OO 1*1 
J q(l-t,x -y)((t,y)dtdy 

/OO r 1 
y g(l - t, x - y)(l - h(y))((t, y)dt dy - p(x) 

with 

/OO /"l 
y q(l-t,x-y)h(y)((t,y)dtdy. 

The measure (1 — h(y))C(t, y)dy is C£°, uniformly in £, so its convolution with the law 
of Ai_( is also smooth, and we obtain that p + p is smooth. It is therefore sufficient 
to study the regularity of p. 

The function g(l —t, .) is smooth out of Z_, so we can differentiate p on (—1/3, 0) 
and obtain 

i.5/3 /.1 

Dp{x)=l I Dq(l-t,x-y)h(y)((t,y)dtdy 

5/3 /-l 

/ (1 - t)e*- 1 J Dg (l - t, x + 1 - y)h(y)C(t, y)dt dy + 0(1) 

as x t 0, because the rest involves the function gm of (1361) which is smooth on 
(—2, 0). The self-similarity of A enables to write 

Dp(x) = J V3 J\l- tf-^e^Dq, (l, %)C(*, y)* ^ + 0(1). 

The law of the stable variable is symmetric and unimodal, so Dqo is nonnegative on 
R_, and is bounded below by a positive constant on some [— C2, — Ci] C (—00, — 2 l l a }. 
Thus the double integral can be bounded below by considering only the part where 
the fraction in Dq is in [— C2, — Ci], where 1/2 < t < 1, and where 1 < y < 4/3. 
With the change s = 1 — t, we obtain 

Dp(x)>c / s^^CCl-s^Jdsdy-C (58) 

with 

J{y) = {0 < 8 < 1/2; d < s-^ a (y -l- x )<C 2 y 
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On the other hand, from (JST1) . 

C(f, y) = E[(y - X,)- 1 -"] > e~*E[(y - X, )- 1 ^] 

where e _i is the probability for A 1 to be up to time t, and X° is the process A 
killed at the exit of M. We have from [7] that the law of X t ° with initial condition 
Xq = is bounded below and above by some positive constants times (1 — u) a ^ 2 du 
for < u < 1 and 1/2 < t < 1, so 

C(l -s,y)>c [ (y- u)- a -\l - uY' 2 du 



> c / ((y - 1) + (1 - «))- a - x (l - M ) a/2 c/M > c'(y - 1)~ Q/2 (59) 

■/l-3( V -l) 

for 1 < y < 4/3 and < s < 1/2. We also have 

J(y) = {s>0;C 1 < s^ a (y - 1 - x) < C 2 ] 

for 1 < y < 5/3 and -1/3 < x < 0, because < y-l-x < 1 and C 2 a < C{ a < 1/2; 
thus 

f s^^ds = C(y - 1 - xf a ~ 2 . (60) 
It follows from (15^,(159"]) and p)]l that 

M/3 

Dp(x) > cy (y - 1 - x) 2a - 2 (y - l)- a/2 dy - C. 

If a < 2/3, we obtain that Dp(x) tends to +oo as x | 0, so Dp(x) tends to — oo and 
the law of Xi is not C 1 . 
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